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Abstract 

From the topological properties of a three dimensional vector order parame- 
ter, the topological current of point defects is obtained. One shows that the 



> 

■^j- ■ charge of point defects is determined by Hopf indices and Brouwer degrees. 

The evolution of point defects is also studied. One concludes that there exist 

o 



crucial cases of branch processes in the evolution of point defects when the 
Jacobian D(&) = 0. 

PACS number(s): 03.40.Kf, 11.27.+d, 47.32. Cc, 02.40.Pc 

I. INTRODUCTION 

In recent years, an outstanding development in the theory of condensed matter is to 
study the defect by topology. It has provided new insights and spectacular predictions. In 
particular there has been progress on the study the defects associated with a nonconserved 
n-component vector order parameter field (f>(f,t) |j],0. In studying such objects in physics 
arise to how one can define quantity like the density of defects. Here comes the important 
conjecture proposed by Halperin, Liu and Mazenko [Hffl, the density of such system has 
been wrote as 
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p = 5($)D( ( f ) /x), (1) 

which is the fundamental equation for such problem. 

In this paper, we will present a new topological current involved in 4-dimensional system 
by use of the 0-mapping topological current theory, which is important in studying the topo- 
logical invariant and structure of physics systems and has been used to study topological 
current of magnetic monopole [|], topological string theory ||, topological characteristics 
of dislocations and disclinations continuum J7|, topological structure of the defects of space- 
time in the early universe as well as its topological bifurcation [|]||, topological structure of 
Gauss-Bonnet-Chern theorem [|3~0| , [3~TH and topological structure of the London equation in 



superconductor JT^] . It is shown that the topological charges of defects can be classified by 
Brower degree and Hopf index of the 0-mapping. It must be pointed out that the existence of 
this topological current is inevitable and necessary, which carries all the important topologi- 
cal properties of this physical system, and includes the defect density even other conjectures 
given by Liu and Mazenko [|J . The further research shows that there exist the crucial cases 
of branch process in the evolution of the point defects when the Jacobian D(&) = 0. We 
calculate the different branches of the worldlines of defects by using the implicit function 
theorem. It is found that the worldlines of point defects will split or merge at the critical 
points. 

II. THE TOPOLOGICAL CURRENT OF POINT DEFECTS 

Consider a 4-dimensional system with a nonconserved 3-component vector order param- 
eter field cj)(r,t). In order to compose a topological current, we introduce a unit vector 
field 

d> a 
n a = JL- y a = 1,2, 3 (2) 

where 
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(x) is just the vector order parameter. The topological current of this system is given by 



AV — (^HV 1 v 2 v 3( - o n ai r) n a2 r) n a3 

J „ ^ ^ai 02013 '~'v\"> '-'V2 w i>z ,c 

on 



It is clear that the topological current is identically conserved, i.e. 



(3) 



d,f = 0. 



(4) 



The topological charge density correspondingly will be defined as 



P = 3 



(5) 



So, Eq. (f|) is just the continuity equation of this system. Considering the following formula 



<9X = \ d ^ a + <^T> 



we can write Eq. (J3|) as 



1 d d 1 

7 ^= F W1M2M3 F f) Affi A0.2X (k a 3_ _C ^ 



in 



If we define 



-a±a 2 a 3 p>/i 
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ze ^ M a^d^d^, 



in which J°(-) is just the usual 3-dimensional Jacobian determinant 



D o ( r ) = D rL) 



X 



X 



I dtp 1 dtp 1 d^ 1 \ 

dx 1 dx 2 dx 3 

d<f> 2 d<t> 2 d<t> 2 

dx 1 dx 2 dx 3 

d<t> 3 d<j> 3 d<f> 3 

\ dx 1 dx 2 dx 3 I 

and make use of the 3-dimensional Laplacian Green's function relation 
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A, 



-An5\<p) 



where A^ is the 3-dimensional Laplacian operator in <p space, we do obtain the 5 function 
like current 



f = 5 3 ($)D^). (6) 

,i 

When \i = 0, we know that the density of j M can be written as 

p = 5\$)D{^). (7) 

From this expression, we find that J M does not vanish only when = 0, i.e. 

<p a (t,x 1 ,x 2 ,x 3 ) = 0, a = 1,2, 3. (8) 

Suppose that the vector field possesses I zeroes (i.e., there exist / point defects) denoted 
as Zi(i = 1, ..., I ). According to the implicit function theorem [I3|, when the zeros points z; t 
are the regular points of <fi which requires the Jacobian 



^©U = ^°©U^o- (9) 

Jb Jb 

the solutions of Eq. (§) can be generally obtained: 

x = Zi{t) i = 1,2, ...,/ (10) 

which represent the worldlines of I point defects moving in space. From Eq. (||), it is easy 
to prove that 

D*(*)\,. = D{*)\ M .*f. (11) 

V s V % dt y ' 

So, the topological current (|3|) can be rigorously expressed as 

j"=mn(^ (12) 

As we proved in Ref. ||, we have 

^ = £/W 3 (2-*f|-L, (is) 

i=i at 

where the positive integer $ is the Hopf index and rji = signD(<fr/x) Zi = ±1 is Brouwer 
degree [|,0. According to Eq. (|T3|), the density of topological charge can be rewritten as 



I 
p = J o = Y / P t r ll S 3 (x-z 1 ). (14) 

It is clearly that the inner structure of this system is characterized by Hopf indices $ and 
Brouwer degrees rji, which are topological invariants. From discussions above, we see that the 
density p(x) is similar to a system of I classical point-like particles with topological charge 
PiTji moving in the 4-dimensional space-time. The topological charge $7^ is also called the 
topological charge of the ith point defect and p(x) can be regarded as the density of defects. 
And , solution ( |T0D can be regarded as the trajectory of the zth defect. From formula (|H"), 
we obtain that the total charge of the system 

1 

Q = I p(x)d 3 x = Y, PiVi- ( 15 ) 

i=i 

The result (@) also has been carried out by Halperin ||, and exploited by Liu and 

Mazenko ||: the first ingredient is the rather obvious result 

Y,6(r-f a (t)) = 8($(r,t))\Dfy\ 

a X 

where the second factor on the right-hand side is just the Jacobian of the transformation 
from the variable to r. This is combined with the less obvious result 

7] a = sgnD((f)/x) \? a 
to give 

P(r,t) =Y,r)J(f- r a (t)) = 8{$)D{4>/x). (16) 

a 

Here we see that the result ([TB]) obtained by Halperin, Liu and Mazenko is not complete. 
They did not considering the cases fa 7^ 1 and D((f>/x) = 0, i.e., r\i is indefinite. It is 
interesting to discuss what will happen and what does it correspond to in physics when 
D((j)/x) = 0. 



III. THE BIFURCATION OF WORLDLINES OF POINT DEFECTS 

As being discussed before, the zeros of the smooth vector (locations of defects) play 
important roles in studying the evolution of point defect. In this section, we will study the 
properties of the zero points, in other words, the properties of the solutions of following 
equations 

^(a; ,^ 1 ,^ 2 ,^ 3 ) = 

f(x°,x\x 2 } x 3 )=0 ■ (17) 

j) 3 (x°,x l ,x 2 ,x 3 ) = 

As we knew before, if the Jacobian determinant 

we will have the isolated solutions fllCf ) of Eq. ([TJ). The isolated solutions are called regular 
points. It is easy to see that the results in section II is based on this condition. However, 
when this condition fails, the above results will change in some way, and will lead to the 
branch process of topological density and give rise to the bifurcation. We denote one of zero 
points as (£*, zi). 

A. the branch process at a limit point 

It is well-known that when the Jacobian D( ^)|(t*,z i ) = 0, the usual implicit function 
theorem is no use. But if the Jacobian 



> x /w,*i) d(t,x 2 ,x 3 ) 






we can use the Jacobian D 1 (-)|( t * ^ instead of D{^)\i t *^ : for the purpose of using the 
implicit function theorem [I3|. Then we have an unique solution of Eqs. ( |TTD in the neigh- 
borhood of the point (£*, Zj) 

t = t(x l ) 



x l = x\x l ) z = 2,3. (18) 

with t* = t(x 1 ). And we call the critical points (t*,Zi) the limit points. In the present case, 
it is easy to know that 



^i ^(i;i(r,*) 

dt l{t *' Zt) D 



(**.*) - n/^M = °° ^ 19 ) 

(t*,Zi) 



i.e. 

dt 

dx 1 



(t*A) - °- 



Then we have the Taylor expansion of Eq. (|T^) at the point (t*, z\) 
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x — z 4 



t = * + d£ l <r* ){x ~ Zl) + 2Wf kt ' A) 

1 d 2 t 
= t * + 2WT lit * A){x ~ Zi) ■ 

Therefore 

which is a parabola in the x 1 — t plane. From Eq. (|20|), we can obtain two worldlines of 
two point defects x\(t) and x^t), which give the branch solutions of the system (p!7|). If 
t^tv2 \(t*,zi) > 0, we have the branch solutions for t > t*, otherwise, we have the branch 
solutions for t > t* . It is clear that these two cases are related to the origin and annihilation 
of point defects. Since the topological charge of point defect is identically conserved, the 
topological charge of these two must be opposite at the zero point, i.e. 

PaVa = -pi2Vi2- (21) 

One of the result of Eq. (|T9|), that the velocity of point defects is infinite when they are 
annihilating, agrees with that obtained by Bray [14]] who has a scaling argument associated 
with point defects final annihilation which leads to large velocity tail. From Eq. ([TT?D, we 
also get a new result that the velocity of point defects is infinite when they are generating, 
which is gained only from the topology of the three dimensional vector order parameter. 



For a limit point, it also requires the D 1 ^-)]^*,^) 7^ 0. As to a bifurcation point ||15|| , 
it must satisfy a more complement condition. This case will be discussed in the following 
subsections in detail. 



B. the branch process at a bifurcation point 



In this subsection, we have the restrictions of the system fll7|) at the bifurcation point 



(t*,zi): 



0(f) It**) 











(22) 



These will lead to an important fact that the function relationship between t and x 1 
is not unique in the neighborhood of the bifurcation point (zi,t*). It is easy to see from 
equation 



dx 1 
~dt 



D 1 



\{t*,z l y- 



\(t*,Zi) 



D{*)k*A) 



(23) 



which under the restraint (p2|) directly shows that the direction of the worldlines of point 
defects is indefinite at the point (z\,t*). This is why the very point (zi,t*) is called a 
bifurcation point of the system ( |17) . 

Next, we will find a simple way to search for the different directions of all branch curves 
at the bifurcation point. Assume that the bifurcation point (z\,t*) has been found from 
Eqs. ([H]) and (|22|) . We know that, at the bifurcation point (zi, £*), the rank of the Jacobian 



matrix [#] is smaller than 3. First, we suppose the rank of the Jacobian matrix [^] is 2 
(the case of a more smaller rank will be discussed later). Suppose that the 2x2 submatrix 



Ji( 



is 



X 



I dtp 1 dcj) 1 \ 

dx 2 dx^ 

a<t> 2 a<f> 2 

\ dx 2 dx 3 ) 



(24) 



and its determinant D l (-) does not vanish. The implicit function theorem says that there 
exist one and only one function relation 
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x* = f(x\*), z = 2,3 (25) 

We denoted the partial derivatives as 

# = ?£. P = d _L. P ° 2 f l . P ° 2 f l . ^ d f* 

From Eqs. (pj]) and (p5|) we have for a = 1, 2, 3 

(f> a = P(x 1 ,f(x 1 ,t),f(x 1 ,t),t) = (26) 



which give 



<;+ES£i£ = ° <») 



dx 1 ^ <9p 9a: 1 



J=2 

^° - ^ + V dr dfj - o ran 

from which we can get the first order derivatives of f l : f{ and f\. Differentiating Eq. ([27D 
with respect to x 1 and £ respectively we get 

e ^°/ii = - jimji + EWkfbfii - #i « = !> 2 ' 3 ( 29 ) 

3=2 3=2 fc=2 

E #/£ = -izW jt f{ + <%ifl + J2(4>%ft)f(] ~4>u a = 1, 2, 3 (30) 

3=2 3=2 k=2 

And the differentiation of Eq. ( pHj ) with respect to t gives 

E #/£ = - Et 2 ^/' + E %/*)/* 1 - ^ a = x > 2 ' 3 ( 31 ) 

3=2 3=2 fc=2 

where 



(32) 



^' fc <9x^x fc ' ^ dxWt 

The differentiation of Eq. ( pS| ) with respect to x 1 gives the same expression as Eq. (|30|). By 
making use of the Gaussian elimination method to Eqs. (§^), Q3"0l) and ([H]) we can find the 
second order derivatives f\ x , j\ t and j l tt . The above discussion does no matter to the last 



component 3 (x, t). In order to find the different values of dx 1 jdt at the bifurcation point 
(zi,t*), let us investigate the Taylor expansion of 4> 3 (x,t) in the neighborhood of (zi, t*). 
Substituting Eq. ( pop into <p 3 (x, t) we have the function of two variables x 1 and t 



F(x 1 ,t)=^(x 1 ,f 2 (x\t),f(x 1 ,t),t) 



(33) 



which according to Eq. ([17]) must vanish at the bifurcation point 



F(zin = o. 



(34) 



From Eq. ( |33"D we have the first order partial derivatives of F(pc l ,t) 
dF 



dx 1 



j=2 Ul j=2 



Using Eqs. 



and (p8|) the first equation of fl22|) is expressed as 

3 



D(Z) 



x 



m,t*) 



E <f>)fi 

3=2 


(j>\ 


<f>l 


E <P)f{ 
3=2 


<f>l 


<& 


<P\ 


03 


01 







(Zi,t*) 



which by Cramer's rule can be written as 



Dt 



dF 



x dx 



r detJi(-)| (W)t .) = 



Since det Ji(|)|(z^*) ^ 0, the above equation gives 

dF , 



With the same reasons, we have 



dx 1 



dF 
~dt 



\(zi,t*y- 



o. 



\{z l ,t*y- 



(35) 



(36) 



(37) 



(38) 



The second order partial derivatives of the function F are easily to find out to be 



d 2 F 

{dx 1 )'- 



3=2 k=2 



(39) 
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SL = k + h<t>y? + 4/1 + #/i + h<p%ft)fi] (40) 

UX Ul j=2 k=2 



— ■ ^ + iiMtfi + flf* + Eitikfftfi] (4i) 



UL j=2 fc=2 

which at (^,t*) are denoted by 

d 2 F d 2 F d 2 F 

A = W? l( * ,t-) ' B = d^dt '<*•*•>' ° = ~W '«'*•) • (42) 

Then taking notice of Eqs. ( |34j), (|37|) , ( p8|) and ([|I|) the Taylor expansion of F(x x , t) in the 
neighborhood of the bifurcation point (zi, t*) can be expressed as 

F{x\t) = \A{x x - z\f + Bix 1 - z]){t - f) + \c(t - t*f (43) 

which by Eq. Q3"3] ) is the expression of <fr 3 (x, t) in the neighborhood of (zi, t*). The expression 
(filf) is reasonable, which shows that at the bifurcation point (zi,t*) one of Eqs. (|l7|). 
4> 3 (x,t) = 0, is satisfied, i.e. 

A{x l - z]f + 2B(x 1 - z]){t - f) + C(t - t*) 2 = 0. (44) 

Dividing Eq. (|44]) by (t — t*) 2 and taking the limit t — > t* as well as x l — > z\ respectively 
we get 

^r> 2 + u 4 + ^°- < 45 > 

In the same way we have 

where A, B and C are three constants. The solutions of Eq. (|4"5| ) or Eq. ( [46]) give different 
directions of the branch curves (worldlines of point defects) at the bifurcation point. There 
are four possible cases, which will show the physical meanings of the bifurcation points. 

Case 1 (A^O): For A = A(B 2 - AC) > from Eq. (g|) we get two different directions 
of the velocity field of point defects 

11 



dx 1 , -B ± s/B 2 - AC , , 

■* L2= 1 ' (47) 

where two worldlines of two point defects intersect with different directions at the bifurcation 
point. This shows that two point defects encounter and then depart at the bifurcation point. 



Case 2 (A ^ 0): For A = 4(B 2 — AC) = from Eq. (45) we get only one direction of 
the velocity field of point defects 

dx 1 . B , 4 . 

or '«= -x (48) 

which includes three important cases, (a) Two worldlines tangentially contact, i.e. two 
point defects tangentially encounter at the bifurcation point, (b) Two worldlines merge into 
one worldline, i.e. two point defects merge into one point defect at the bifurcation point, (c) 
One worldline resolves into two worldlines, i.e. one point defect splits into two point defects 
at the bifurcation point. 

Case 3 (A = 0, C ^ 0): For A = 4(5 2 - AC) = from Eq. (g§) we have 



dt . -B ± V-B 2 - AC n 2B 

3?l"= c = °- ~c- (49) 

There are two important cases: (a) One worldline resolves into three worldlines, i.e. one 
point defect splits into three point defects at the bifurcation point, (b) Three worldlines 
merge into one worldline, i.e. three point defects merge into one point defect at the bifur- 
cation point. 

Case 4 (A — C — 0): Eq. fl45|) and Eq. fl46|) gives respectively 



— = 0, — = 0. 50 

dt dx 1 y J 

This case is obvious similar to case 3. 

The above solutions reveal the evolution of the point defects. Besides the encountering 
of the point defects, i.e. two point defects encounter and then depart at the bifurcation 
point along different branch curves, it also includes spliting and merging of point defects. 
When a multicharged point defect moves through the bifurcation point, it may split into 
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several point defects along different branch curves. On the contrary, several point defects 
can merge into one point defect at the bifurcation point. The identical conversation of the 
topological charge shows the sum of the topological charge of final point defect (s) must be 
equal to that of the initial point defect(s) at the bifurcation point, i.e., 

EM/ = Eta ( 51 ) 

/ 

for fixed I. Furthermore, from above studies, we see that the generation, annihilation and 
bifurcation of point defects are not gradual changes, but start at a critical value of arguments, 
i.e. a sudden change. 

C. the branch process at a higher degenerated point 

In above subsection, we have studied the case that the rank of the Jacobian matrix Nj2] 
of Eqs. fll7f) is 2 = 3 — 1. In this subsection, we consider the case that the rank of the 
Jacobian matrix is 1 = 3 — 2. Let the J 2 { ) = afr anc ^ suppose that det J 2 7^ 0. With the 
same reasons of obtaining Eq. fl25|) , we can have the function relations 

x* = f\x\x\t) . (52) 

Substituting relations ( |52] ) into the last two equations of (|T7D, we have following two equa- 
tions with three arguments x 1 , x 2 , t 

F^x^x 2 ^) = <f) 2 {x 1 ,x 2 ,f 3 {x 1 ,x 2 ,t),t) = 

(53) 

F 2 {x 1 , x 2 , t) = ^{x 1 , x 2 , / 3 (x\ x 2 , t),t) = 0. 

Calculating the partial derivatives of the function iq and F2 with respect to x 1 , x 2 and t, 
taking notice of Eq. ( p2[ ) and using six similar expressions to Eqs. fl3~T|) and ([38|), i.e. 

aF * l - n 9F J l - n dF * l - n ,--12 f*rt 

we have the following forms of Taylor expressions of F% and F 2 in the neighborhood of (,?;, t*) 

F {x\ x 2 , t) « A^x 1 - z\f + A^x 1 - z]){x 2 - z 2 ) + A^x 1 - z\) 
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*\2 



x(t - f) + A j4 (x 2 - zfY + A^(x 2 - zf){t - f) + A j6 (t - t*) 



J = 1,2. 



(55) 



In case of A j± ^ and A,- 4 ^ 0, dividing Eq. fl55|) by (t — t*) 2 and taking the limit t — > t*, 



we obtain two quadratic equations of %- and dl " 



dt 



Ax 1 
~dt 



1 Jrr2 



dx dx 






Aji ^_y + Aj2 __ + A . 3 _ + A u (—) 2 H .4 ,,-^ ! .4 /(1 



dx 2 



(56) 



J = 1,2. 



Eliminating the variable dx 1 /dt, we obtain a equation of dx 2 /dt in the form of a determinant 



A n A 12 v + A 23 A l4 v 2 + A 15 v + A w 

A 1X A 12 v + A 13 A 14 v 2 + A 15 v + A 16 

A 21 A 22 v + A 23 A 24 v 2 + A 25 v + A 26 

A 21 A 22 v + A 23 A 24 v 2 + A 25 v + A 



26 



(57) 



where v = dx 2 /dt, that is a ^th order equation of dx 2 /dt 

ao(-i-) 4 + «i(-ir) 3 + fl 2(--) 2 + a 3 — - + a 4 = 0. 
at at at at 



(58) 



Therefore we get different directions of the worldlines of point defects at the higher degen- 
erated point bifurcation point. The number of different branch curves is at most four. 

At the end of this section, we conclude that in our theory of point defects there exist the 
crucial case of branch process |lTJ. Besides the encountering of the point defects, i.e. two 
point defects encounter and then depart at the bifurcation point along different worldlines, 
it also includes spliting and merging of point defects. When a multicharged point defect 
moves through the bifurcation point, it may split into several point defects along different 
worldlines. On the contrary, several point defects can merge into one point defect at the 
bifurcation point. Since the topological charges of point defects is identically conserved (|]), 
the sum of the topological charge of final point defect (s) must be equal to that of the initial 
point defect (s) at the bifurcation point. 
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IV. CONCLUSIONS 

We have studied the evolution of the point defects of a three dimensional vector order 
parameter by making use of the mapping topological current theory. We conclude that 
there exist crucial cases of branch processes in the evolution of the point defects. This means 
that the point defects generate or annihilate at the limit points and encounter, split or merge 
at the bifurcation points of the three dimensional vector order parameter, which shows that 
the point defects system is unstable at these branch points. There are two restrictions of the 
evolution of point defects in this paper. One restriction is the conservation of the topological 
charge of the point defects during the branch process, the other restriction is the number 
of different directions of the worldlines of point defects is at most four at the bifurcation 
points. At last, we would like to point out that all the results in this paper are obtained 
from the viewpoint of topology without using any particular models or hypothesis. 
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